HEAVY FLAVOR PHYSICS : LIFETIMES AND FLAVOR CHANGING 
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The Tevatron, with pp collisions at yfs =1.96 TeV, can produce all flavors of B hadrons and 
allows for unprecedented studies in the B physics sector. These range from measurements of 
B hadron properties to searches of new physics signatures. The CDF'— ' and DO'-' detectors 
currently have more than 7 fb _1 of data recorded. This paper presents a selection of recent 
results on lifetimes and flavor changing neutral currents using between 4.3-5.0 fb _1 of data. 



1 B hadron lifetimes 

The experimental measurement of B hadron lifetime ratios is an important test of the theoretical 
approach to B hadron observables known as the heavy quark expansion. The ratio t(B + )/t(B°) 
(charge conjugates are implied throughout) is predi cted^ ' ^ I to be in the range 1.04-1.08 and 
the ratio t(A^)/t(B°) in the range 0.83-0.95. ' The measured world average B + and B° 
lifetimes are dominated by the Belle experiment^ Of recent interest is the A° lifetime. Until 2006 
all measurements were in agreement but lay at the lower end of the theoretically expected v alue. 
Since then, two high precision CDF measurements are significantly above previous result s. ^ I 
The analysis described here is the most precise measurement of the B + , B°, and A° lifetimes 
/^ • and ratios. 

The B + , B°, and A& lifetimes were measured using 4.3 fb _1 of data with decay channels 
B + -> J/ipK + , B° — > J/ipK* , B° — > J/ipK s and A& — > J/ipA. In previous measurements the 
uncertainty due to detector resolution has been a leading source of systematic uncertainty. In 
this analysis the proper decay time is determined using the J/tp vertex to provide similarity in 
the decay time resolution between channels and to allow for the cancellation of certain systematic 
uncertainties. A detailed resolution model is also introduced in this analysis. The signal decay 
time is modelled as an exponential decay convolved with the resolution model. The resolution 
model is a superposition of three Gaussians. They are each centred at t = 0, and have a 
width of event decay time uncertainty, of 4 , multiplied by a scale factor. The restriction to 
models symmetric about t = is motivated by simulation, while the number of components is 
determined from data. The parameters of the resolution function are determined from the mass 
sidebands as the fraction of background events expected to originate from the primary vertex is 
between 80-90%, depending on channel and background model, and therefore provides a useful 
sample from which to determine the resolution. The overall fit is an unbinned likelihood fit to 
the mass, decay time and decay time uncertainty distributions simultaneously. The projections 
of the mass and decay time distributions from the A& data are shown in Fig. [TJ 

We measure t b + = 1.639 ± 0.009 (stat) ± 0.009 (syst) ps, t b o = 1.507 ± 0.010 (stat) ± 
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Figure 1: The A& mass and proper decay time distributions. 
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1.537 ± 0.045 (stat) ± 0.014 (syst) ps. The lifetime ratios are cal- 
culated as t b+ /t b o = 1.088 ± 0.009 (stat) ± 0.004 (syst) and t a o/t b o = 1.020 ± 0.030 (stat) ± 

. . b 

0.008 (syst).'^ These are the world's best measurements of the lifetimes and ratios. The im- 
provement in the systematic uncertainty from 0.033 ps (1.0 fb -1 to 0.014 ps (4.3 fb -1 ) is evident 
in the t(A&) measurement. The A& lifetime remains higher than the world average but is not 
inconsistent with theoretical predictions. 

Measurements of other B baryon lifetimes are also interesting as they can va lidate the 

-I O 1 o 

predi ction o f the baryon hierarchy and ratios. CDF and DO have both observed the Hj , * ' and 
the Q b ' ' baryons. CDF has performed lifetime measurements on these samples. The low 
statistics, particularly in background motivates a different approach to lifetime measurement. 
The data are divided into bins of proper decay time. A mass fit in each bin determines the 
number of signal candidates in each bin. The lifetime is then determined from the yield in each 
bin compared to the expected distribution for a given lifetime. The method has been validated 
on Afe and B° decay modes and gives results consistent with those from other methods. Using 
4.3 fb" 1 CDF measures t(E{,) = 1.56^Q 2 5(stat) ± 0.02(syst) ps which is the first measurement 
using a fully reconstructed decay of E&, and r(ilfe) = 1.13l ' 40 (stat) ±0.02(syst) ps which is the 
first measurement of r(Q&jr"H 



2 Flavor Changing Neutral Currents 

Flavor changing neutral currents (FCNC) are forbidden at tree-level in the standard model. 
Processes that involve FCNC through higher order diagrams can provide information on new 
physics. New particles participating in the higher order diagrams can enhance decay rates or 
alter expected kinematic distributions from the standard model predictions. FCNC provide a 
complementary approach to new physics searches alongside direct searches, as any signatures 
of new physics that are observed give information on the flavour structure of new physics, and 
thus constrain the set of new physics models that are consistent with data observations. 



2.1 B 
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One highly studied decay is B s — > fi + fi~. The standard model branching fraction is^ (3.42 ± 
0.54) xl0~ 9 which is beyond the CDF and DO detector experimental sensitivity.However, these 
branching ratios can be enhanced by a factor of 10 — 100 by super symmetric models or other new 
physics. The decay B° — > /^ + /i~ is further suppressed by the ratio of CKM elements, |14rf/T4 s | 2 . 
Given the detector sensitivities, observation of either decay is an unequivocal signature of new 
physics. If observed, the ratio of measured B s and B° branching fractions would give information 
on the flavor structure of the new physics. 



The analyses carried out at CDF and DO are similar; the CDF analysis uses 3.7 fb _1 and is 
described below. One challenge in this analysis is to reduce the large backgrounds. A number 
of baseline selection requirements are applied which result in a reduction of background by a 
factor of 300 while 50% of signal would remain in the geometric and kinematic acceptance of the 
detector. The baseline selection includes selection on transverse momentum, vertex quality, and 
muon ID algorithms which reduce backgrounds where hadrons have been misidentified as muons 
and decays involving a kaon that has decayed inflight to a muon. For further enhancement of 
signal events a neural network is used. It is trained on six variables; the proper decay time 
and proper decay time significance, the transverse momentum of the di-muon candidate, the IB- 
candidate track isolation, the px of the lower momentum muon candidate, and the 3D opening 
angle between the vectors yp^ and the displacement vector between the primary vertex and the 
dimuon vertex. The neural network is trained using signal events generated by Monte Carlo 
and mass sideband events for background. The neural network output is shown in Fig. [2j 
The remaining background is estimated from continuum combinatorics from sidebands and mis- 
reconstructed B — > hh decays which peaks in the signal region. The background estimates 
are cross-checked using control samples from data such as like-sign muons. The background 
predictions are compared with the observed data and no statistically significant discrepancies 
are observed. 

A relative normalisation to the channel B + — > J/ipK + is used to determine the B — > fj,fi 
branching fraction. The number of observed signal events, N s , and normalization events, N + , 
can be used to obtain the branching fraction via: 

B(B° a -+ ,iV~) = #■•-• T ■ ®( B+ "> J/4>K + , J/V -)• /xV), (1) 

JV + € s Js 

where f u , f s are fragmentation fractions and e + and e s are the efficiencies for the normalisation 
and signal channels. The efficiencies are determined through a combination of Monte Carlo 
simulation and data driven techniques involving samples of J/ip, B — > J /ijjK + and B — > J/ipcp. 
After selection, the data observed in the mass ranges corresponding to B s and B° are 
consistent with there being only background events and are shown in Fig [2] in three sepa- 
rate bins of the neural network discriminant. CDF sets limits on the branching fractions for 
these decay processes. Using a data sample of 3.7 fb _1 CDF extracts a 95% (90%) C.L for 
B(B S -> ft*/*-) < 4.3 x 10- 8 (3.6) x lO" 8 and B(B° -> fi+n') < 7.6 x 10~ 9 (6.0) x 1(T 9 OH These 
are the world's best limits on these branching ratios. At DO an analysis of 2 fb _1 obtained a 
limit of B(B S -)■ fi + n~) < 9.5 x 10~ 8 (7.5) x 10" 8 .Q2]The analysis of 5 fb _1 is on-going with the 
signal region still blinded. The expected upper limit on the branching fraction on these data 
is B(B S — > p + ijl~) < 5.3 x 10 _8 (4.3) x lO -8 .^^ These measurements at the Tevatron have lim- 
ited the parameter space of allowed new physics models as large enhancements of the B — > fj,fi 
branching ratio are inconsistent with the observed data. 

2.2 B ->■ n+ii~h 

One group of interesting decays to study are B — > [ifih where the B is either B + , B° or B s and h 
stands for either K + , K* or </>, respectively. The branching ratios of these decays are 0(1O -6 ) and 
are observable. New physics processes can enhance decay amplitudes. The interference between 
amplitudes due to new physics particles and the standard model processes may be observable in 
effect on branching ratios, polarisation and forward backward asymmetry. Comparison of these 
observables to standard model expectation can indicate whether the underlying dynamics are 
governed by the standard model or other models such as SUSY, or 4th generation quarks. 

The three decay channels are reconstructed from data collected by triggers that require 
two charged particles with px > 1.5GeV/c or 2.0 GeV/c with corresponding hits in the muon 
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Figure 2: The left plot shows the neural network discriminant for signal and background. The right plot shows 
the data after selection in three bins of the neural network discriminant. 



chambers. Vertices comprising of the muon pair and the h are required to have a x 2 probability 
of at least 10 -3 . For the normalisation of the branching ratios the control channels B — > J /i^h 
are also reconstructed. The same final state allows a number of systematic uncertainties to 
cancel. Selection criteria are placed on the decay time significance, transverse momentum, 
impact parameter of the B candidate, the eft and K* mass, and particle identification to reduce 
the combinatoric and fake muon backgrounds. The normalisation channels are required to 
have the di-muon invariant mass within 50 MeV/c 2 of the J/ip mass. For the signal channels, 
candidates with a di-muon mass in the region near the J/ip and tp' are rejected. Other criteria 
are applied to reduce other peaking backgrounds. After the loose selection a neural network 
technique is used to provide the final selection. It is optimised to maximise both the branching 
ratio and the lepton forward-backward asymmetry, ApBi significance. 

The signal yield is determined from an unbinned maximum likelihood fit to the B invariant 
mass distribution. The signal shape is two Gaussians with different means, and is determined 
from Monte Carlo with the B mass resolution scaled to data in the normalisation channels. The 
background is a first or second order polynomial. The mass distribution and yields are shown in 
Figure [3l For each channel a significance greater than 6<r is observed. The measured branching 
ratios are B{B + -r K + fi + n~) = [0.38 ± 0.05(stat) ± 0.03(syst)] x 10" 6 , B(B° -> K * /i+//~) = 
[1.06±0.14(stat)±0.09(syst)]xl0- 6 , B{B S -)• <f>(i + fT) = [1.44±0.33(stat)±0.46(syst)] x lO" 6 .^ 
These numbers are consistent with previous results and other B-factory measurements. In the 
case of B s — > (frfi/i, this is the first observation of this decay channel and is also the rarest B s 
decay observed so far. 

The differential branching ratio as a function of the di-muon invariant mass is determined 
by binning the data into 6 bins and repeating the branching ratio measurement. The 6 bins 
are chosen to correspond to the analysis performed at Belle to aid comparison between the two 
experiments. The results are shown in Fig. [4) The region consistent with the standard model 
is between the two lines; this band exists due to uncertainty in form factors. No inconsistency 
with the standard model is observed. 

The B + and B° decays are used to determine Afb and the longitudinal polarisation of the 
K*. The theoretical uncertainty on the standard model prediction has very small uncertainty 
and hence provides an interesting opportunity to observe the effects of new physics. Aps is 
determined from the distribution of cos(6^), where 6^ is the helicity angle between the n + (fi~) 
direction and the opposite of the B{B) direction in the di-muon rest frame. The polarisation 
is measured from Ok where this is the angle between the kaon direction and the direction 
opposite to the B meson in the K* rest frame. The fit results are shown in Fig. EJ where one 
possible beyond the standard model scenario is also shown. The results are consistent with the 
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Figure 3: The mass distribution of B — >• fifih candidates after selection. 
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Figure 4: The differential branching ratios for B + (right) and B°(left). The hatched regions are charmonium veto 

regions. Solid lines are the SM expectation. 
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Figure 5: Results for Afb and K* polarisation. The solid (dotted) line is the standard model (a new physics 

model) expectation. 

standard model prediction^ and are consistent and of similar precision to those observed at the 
B factories. 

3 Conclusion 

The Tevatron continues to make significant contributions to the field of flavor physics. The 
results presented here have shown the possibilities of hadron colliders for measurements in flavor 
physics despite the high combinatoric backgrounds. The contribution to lifetimes will be used to 
better understand the interaction of quarks inside hadrons. The indirect searches of new physics 
through flavor physics observables remain complementary to direct searches. 

The Tevatron accelerator continues to provide a high luminosity and these analyses can be 
expected to be updated with 10 fb _1 . This will allow for significantly improved measurements for 
the measurements that are currently limited by statistical uncertainty. This will be particularly 
interesting for B — >• fj.fi and B — > fifth where hints of new physics could be observed. 
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